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We study the effect of mode-coupling on a single-photon device in which a dipole-quantum-emitter (DQE) is 
embedded in a bimodal whispering-gallery-mode cavity (WGMC). A scatterer is used to induce mode coupling 
between counter-clockwise and clockwise propagating light fields, which interact with the DQE. In contrast 
to models for the interaction between a DQE and a (one-mode or two-mode) cavity field, we find that strong 
photon antibunching can occur even for a weak DQE-field coupling and large dephasing of the DQE, when 
mode coupling is introduced. We also find that mode coupling can make the device robust against either the 
frequency mismatch between cavity modes and the DQE or the coupling strength mismatch between the DQE 
and each mode in the two-mode cavity. Moreover, we find that these mismatches can be used to generate better 
antibunching in the weak DQE-field coupling regime. Our study shows that mode coupling in a bimodal cavity 
is very important for the realization of a good single-photon device. 

PACS numbers: 42.50.Ct, 42.50.Ar 


I. INTRODUCTION 

Single-photon devices play a crucial role in quantum infor¬ 
mation science, including quantum communication and quan¬ 
tum computing technology ill |2]. These devices are useful 
for quantum key distribution iS® , generation of entangled 
states iH , quantum metrology BH], single-photon quantum 
memory |[1Q| - [1^ . linear optical quantum computing 1114 fl^. 
and quantum simulations |[l^ Ol- Thus, a highly-efficient 
method for generating single-photons with low error is ur¬ 
gently required for quantum information processing. Simi¬ 
lar to Coulomb blockade for electrons in mesoscopic elec¬ 
tronic devices photon blockade, where the subse¬ 

quent photons are prevented from resonantly entering a cavity, 
is a promising way to produce single-photon sources. Photon 
blockade requires that the strength of the photon-photon Kerr 
nonlinear interaction imm is much larger than the decay 
rate of the cavity field. 

In principle, the photon blockade can be measured by the 
second-order correlation functions of photons Pho¬ 

ton blockade in the optical frequency domain has been ob¬ 
served in: a trapped atom coupled to a cavity field 
a two-level system coupled to a bimodal microtoroidal res¬ 
onator ii, and a quantum dot coupled to single-mode field of 
a photonic crystal resonator . Recently, photon blockade 
in the microwave frequency domain was also observed 
by coupling a superconducting artificial atom llsol - ls^] to a 
transmission line resonator. A crucial prerequisite for photon 
blockade in these experiments ll25l - [^ is to reach the strong 
field-atom coupling regime ll^ . i.e., the coupling strength be¬ 
tween the cavity field and the dipole-quantum-emitter (DQE) 
exceeding the decay rates of both the cavity field and DQE. 
The strong DQE-field coupling can induce a strong nonlinear 
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photon-photon interaction and could be used to produce pho¬ 
ton blockade and single-photon sources ll34 - [^ . However, 
when a cavity field works at the single-photon level, it is not 
easy to achieve a cavity-atom-induced strong photon-photon 
interaction. 

Recently, photonic “molecule” systems, consisting of two 
coupled cavities with Kerr nonlinearities, were proposed to 
generate antibunching photon and photon blockade m. In 
contrast to previous studies ll25l - l2^ . it was found that the pho- 
ton blockade can occur even in a weak field-atom coupling 
regime, when an additional cavity is coupled iS. This new 
mechanism corresponds to a destructive quantum interference 
effect in nonlinear photonic molecule systems . Motivated 
by several studies 11^ 1^. various systems were proposed to 
achieve photon blockade and antibunching photons, such as 
coupled optomechanical systems ||4T| - I^ . dipolariton systems 
formed by cavity fields and excitonsJH, coherent-feedback 
controlled optomechanical systems 11^ . and coupled micro¬ 
cavities with second- or third-order nonlinearities 1^ . 
However, the role of the coupled cavity was not explored in 
detail. Thus it is interesting to further explore this issue as a 
function of the coupling strength between the two cavities. 

Motivated by an experiment on a photon turnstile, in 
which a DQE is embeded in a bimodal whispering-gallery- 
mode cavity (WGMC) we here mainly study the mode 
coupling effect on the nonclassic al prope rties of photons. Eor 
two strongly-coupled resonators ll39l - [4^ . two modes of two 
resonators are coupled to different nonlinear quantum sys¬ 
tems or only one of the cavity modes is coupled to a non¬ 
linear quantum system. However, here: (i) a single WGMC 
is used to support two degenerate counter-propagating cavity 
modes; (ii) the mode coupling or normal-mode splitting be¬ 
tween these two counter-propagating modes is realized by a 
scattering object ll4^l - [^ : (hi) two cavity modes are coupled 
to the same DQE. Compared to the studies where a two-level 
system is coupled to two modes of a cavity here the 

mode coupling is introduced by a scattering object. 
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FIG. 1: (Color online) Schematic of a WGMC-DQE system which is 
driven via an optical fiber taper waveguide with a cavity-waveguide 
coupling rate Ke- The nanoparticle in the WGMC, causes an effec¬ 
tive coupling at a rate J, between the clockwise and counterclock¬ 
wise propagating modes. Such two WGM modes have a DQE-field 
coupling rate g and intrinisc cavity decay rate K,i. The DQE has a 
spontaneous emission rate 7a and a phase dephasing rate 7 p. 


The remaining part of the paper is organized as follows. 
In Sec. II, the theoretical model is introduced. In Sec. Ill, 
the measurements for the mode-coupling are discussed. In 
Sec. IV, the photon blockade is analytically and numerically 
studied in different parameter regimes via the second-order 
correlation function. In particular, we analyze the photon 
blockade in the weak DQE-field coupling regime. We also 
compare our results with those in coupled-cavity systems 
1^ . In particular, we analyze the physical mechanism of the 
photon blockade in the system studied here. In Sec. V, we 
study the robustness of the system on the dephasing of the 
DQE. In Sec. VI, we study the robustness of the system on ei¬ 
ther the frequency mismatches between the cavity mode and 
the DQE or the coupling strength mismatch between the DQE 
and each cavity mode. Conclusions and perspective discus¬ 
sions are finally presented in Sec. VII. 


II. HAMILTONIAN 

As schematically shown in Eig. [T] we study a system that 
consists of a single whispering-gallery-mode cavity (WGMC) 
and a dipole-quantum-emitter (DQE). As studied in Ref. 
the WGMC supports a counterclockwise (CCW) and a clock¬ 
wise (CW) propagating modes which have the same fre¬ 
quency. These two modes are coupled to each other through 
a nanoparticle (called as a scatterer) The DQE is modeled as 


a two-level system, which can be a quantum-dot, an atom, or 
other systems. Here we do not focus on any particular system, 
but consider a generic D QE. The Hamiltonian of the whole 
system is given by 11^ 1^ 1^ 

H = huj (a^a + b^b) -b hujg,(7z + hJ{a^b + b^a) 

+ h [gaa^cr- -|- gbb^cr- + sa ^+ H.c.] , (1) 

where a and b (a^ and b^) are the annihilation (creation) oper¬ 
ators of CCW and CW propagating modes with frequency cc, 
respectively. cr_ and its conjugate operator cr+ are the ladder 
operators describing a DQE with frequency uja, and J is the 
mode-coupling strength between the CCW and CW modes. 
We assume that the DQE and two propagating modes have 
coupling strengths ga and g^. The driving field is assumed to 
be coupled to the CCW mode with the coupling strength 5. 

The mode coupling in the WGMC is usually introduced 
when a scatterer is present. The scatterer can be considered as 
a nanoparticle, which is used to enhance and control |[^ 
the effective coupling strength J between the CCW and CW 
modes. Let us consider a simple example where the nanopar¬ 
ticle is assumed to be a nanosphere with radius R X, where 
A is the wavelength of the light field. Then the response of 
a particle to an electromagnetic field can be calculated by the 
Clausius-Mossotti relation and the effective coupling strength 
J becomes |[^ 

J=-^afoj, ( 2 ) 

in the electrostatic limit, with 

Here / is the mode function of the WGMC, V is the mode vol¬ 
ume of the WGMC, Up is the refractive index of the particle, 
and the surrounding medium is assumed to be air. 

The scatterer-induced mode splitting has been observed ex¬ 
perimentally, and proposed for a highly sensitive and robust 
platform for detecting nanoscale objects lld^l - l^ . As reported 
in m, a nanoparticle of radius 40 nm can introduce an ap¬ 
parent mode splitting with a mode-coupling strength J > k,, 
where hc is the decay rate of the cavity field, producing a 
strong-mode coupling. Erom Eqs. ^ and ([3]), it is clear that 
the mode-coupling strength J can be enhanced by increasing 
the radius of the nanoparticle. The effective coupling strength 
can also be enhanced by increasing the number of nanopar¬ 
ticles 11^ . improving the quality factor of the microres¬ 
onator [[Sfl], introducing gain medium into the cavity (HI], or 
using Raman-gain-induced loss compensation (H]- Here, we 
will not pay more attention to how the nanoparticle will affect 
the system. Instead, we just assume that the mode-coupling 
strength J is a positive and real number, which can be con¬ 
trolled and set to be much larger than the decay rate n of the 
cavity field. We will mainly focus on the study of how the 
mode-coupling strength J affects the single-photon behavior 
in the system. 

In the rotating reference frame at the driving-field fre¬ 
quency ccd, i-e., with a unitary transformation U = 
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FIG. 2: (Color online) Normalized transmission spectra Pt and re¬ 
flection spectra Pr of the coupled WGMC-mode system. The hori¬ 
zontal axis shows the normalized A = cj — cjd, which is the detun¬ 
ing from the driving frequency Also, r is the decay rate of the 
cavity held. Two different parameter regimes are shown: (i) when 
the optical mode-coupling rate J is much larger than the cavity-fleld 
linewidth k, corresponding to the red-dotted curve for Pt and red- 
dashed curve for Pr, with J = 6 as an example; and also (ii) when 
J is comparable to r, corresponding to the black-solid curve for Pt 
and black-dash-dotted curve for Pr, with J = 0.8 /c as an example. 

exp \—iujd{o)cL (Tz)t \, the Hamiltonian in Eq. ([T]) is 

given as 

iTeff = (o) a + 6^6) + hAg^az + hJ{a^b + b^a) 

+ h [g{a^ + 6 ’*')cr_ + + H.c.] , (4) 


operators associated with the CCW mode, CW mode, and the 
DQE with zero mean-value, i.e., (ain) = (bin) = (c^’in) = 0- 
The mode coupling can be clearly observed in the transmit¬ 
ted (Pt) and reflected (Pr) cavity field power. Let us now cal¬ 
culate these through the outputs of the modes a and 6 , respec¬ 
tively. Using the mean field approximation and from Eqs. ([5]) 
and dll), the equations of motion for the mean value of each 
operator can be given as 



(*A + 1) (a) - iJ {b) - ig (cr_) - ie, 

( 8 ) 


(*A + 1) (b) - iJ (a) - ig (ct.) , 

(9) 

1 

II 

1 

(*A + 1 ) (cr) + ig{{a) + ( 6 )) (cr^). 

(10) 


Using the input-output theory (bl, [b^], the output of each 
mode can be obtained as 


(ttout) - /— p (tt) 5 (11) 

(bout) = v^{b). ( 12 ) 

Here, we are only interested in the mode coupling induced 
by the scatterer, thus we can set = 0. Then the dynamical 
equations of the cavity modes are simplified to 


(a) = - 

(*A + 1) (a) - iJ (b) - is, 

(13) 

{b) = - 

(*A + 1) ( 6 ) - iJ (a). 

(14) 


with the detunings A = cc — ccd and Aa = cca — from 
the driving field ccd to the cavity modes cc and to the DQE 
cda- Below, we first discuss the case for Qg, = g\y = g and 
Aa = A (i.e., uj = ujg) in Secs. Ill, IV, and V. The mismatch 
of the coupling strengths ^a and g\y and the mismatch of the 
frequencies uj and cca will be discussed in Sec. VI. 


III. MEASUREMENT OF THE MODE-COUPLING 
STRENGTH 


Let us first study how the mode-coupling strength J could 
be measured. The dynamics of the whole system is described 
by the quantum Langevin equations 


da 

dt 

dt 

dcT- 

dt 


(^A + ^ ~ ^^b — igcF- — is — (5) 

^iA + ^ ~ ~ ~ ( 6 ) 

+ a-+igaa^+igba^-y/yain- (7) 


Here, the decay rates of the two cavity modes are denoted by 
K, which consists of the intrinsic loss Ki and the external loss 
Kq, i.e., K, = Ki Kq. Eor the DQE, the decay rate 7 consists 
of the spontaneous emission rate 7 ^ and phase dephasing rate 
7 p, i.e., 7 = 7 a + 27 p. Here, ain, 6 in and (Tin are the noise 


By solving Eqs. (fTSl) and ([TO in the steady-state with (a) = 
(b) = 0 , we have 


(o) 

ib) 


A-if 
(iA + f)V 
-J 

(iA+f)^ + J2^’ 


(15) 

(16) 


Combining Eqs. STTh and ([TO with the steady-state solutions, 
the normalized transmission power Pt and refiection power 
Pr from the WGMC are given by 


Pt — — K^iout)!^ 

PR = ^\{bout)f 


i ^ \/^ -^f) 

(iA+f)%J 2 

2 

yJ~dJ 

(iA+f)Vj 2 


(17) 

(18) 


In Eig.[2l the transmission Pt and refiection Pr powers are 
plotted as a function of the detuning parameter A = cc — ccd- 
Two different parameter regimes are discussed: (i) J > n 
and (ii) J ^ n. Eor J > it is obvious that a pair of res¬ 
onances appear and are located at A ± J. Actually, these 
resonances correspond to two different standing light waves, 
which are superpositions (i.e. (a + 6 ) /a /2 and {a — b)/ \/ 2 ) 














4 


of the two counter-propagating cavity modes. When the cou¬ 
pling strength J is much larger than k, a well-resolved mode¬ 
splitting can be clearly observed in the transmission or reflec¬ 
tion spectra. However, the distance between two resonances 
becomes close when J k, and is no longer distinguishable 
when J n. 


IV. PHOTON BLOCKADE IN THE WEAK DQE-FIELD 
COUPLING REGIME 


In this section we will study the photon blockade. We 
mainly study how the quantum behavior of the mode flelds 
in the WGMC varies with the change of the mode-coupling 
strength J between the two modes. Compared with previous 
studies on the strong coupling between a two-level atom and a 
single-mode cavity (e.g., in Refs. ll25l - [^ ) or a bimodal cav¬ 
ity without the mode coupling |[?^ [^. we will here focus on 
the weak coupling between a DQE and the cavity modes with 
mode coupling. 


A. Master equation and second-order correlation functions 


By taking the dissipation into account, we can write the 
master equation |[59l - [^ of the reduced density matrix p for 
whole system as 

dp 1 
dt ih 

+ ^{2bpb^ — b^bp — pb^b) 

-h ^{2apa^ — a^ap — pa^a) 

+ p), (19) 

in the Markov approximation. Here, i^eff is given in Eq. 0 . 
We assume that the cavity flelds and the DQE are in a zero- 
temperature environment, to simplify the calculation. In the 
number states basis |na,nb,i) and i'), the formal 

solution of p in Eq. (O can be given as : 


+ —{2apo) — a^ap — pa^a) 


P{t) — ^ ^ ^ ^ ^ P'na,nb,i\n^/ ,n^/ ,i' {'^a' \ • 

( 20 ) 

where ria {ria ') represents the photon number of the CCW 
mode, rib i^b ') represents the photon number of the CW mode 
and i{i') = — represents the excited and ground state of 

the DQE, respectively. 

We can obtain the steady-state solution pss of the reduced 
density operator p by setting dp/dt = 0 , and then the statis¬ 
tical properties of the driven CCW propagating cavity mode 
in the WGMC can be studied via the normalized equal-time 
second-order correlation function 




(a)a^aa^ 

(ata)^ 


aa) 

[Tr(pssO^’*‘G^)]^ ’ 


( 21 ) 



FIG. 3: (Color online) The second-order correlation functions 
(0) of a CCW propagating intercavity field versus the detuning 
parameter A//^ = (ct; — ctid)//^, for a number of different parame¬ 
ters. These correlations are calculated numerically using the master 
equation and also analytically using the Schroedinger equation in the 
steady-state. The (0) is shown as an orange-rhombus solid curve 
for the interaction between a DQE and a single-mode cavity field; 
The (0) is shown as a black solid curve (numerically calculated) 
and a cyan-dotted line (analytically calculated) for the interaction be¬ 
tween a DQE and two-mode cavity fields with J = 0. The pi^^(O) 
is shown as a red-circular solid curve (numerically calculated) and 
a green-triangle solid curve (analytically calculated) for the interac¬ 
tion between a DQE and two-mode cavity fields with J = 30/c. The 
system parameters for this simulation are: n = 407a, g = 207a, 
£ = 7a, and 7a = 1 MHz. 


In this section, we first study the antibunching effect of the 
CCW mode by setting 7 ^ = 0, in order to consider the same 
environmental effect as in previous studies m. The effect 
of dephasing will be studied in Sec. V. We point out that the 
statistical properties of the CW mode can also be studied in a 
similar way as those for the CCW mode. 


B. Numerical calculations 

We first study how the second-order correlation functions 
Pqqw(O) vary with the detuning A = cc — ccd between the 
driving field and the cavity modes, and also show how the sta¬ 
tistical properties of the cavity field are changed by the mode¬ 
coupling strength J. In Eig. [3l Pccw(^) versus the detun¬ 
ing A is plotted in the DQE-field weak-coupling regime for 
a DQE interacting with either: (i) a single-mode cavity field, 
or (ii) bimodal cavity fields without mode coupling, or (hi) 
bimodal cavity fields with mode coupling. Eigure [3] clearly 
shows that there is only photon bunching when only a single¬ 
mode cavity field is coupled to the DQE. This is due to the 
weak DQE-field coupling, i.e., ^ < n. When the two modes 
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FIG. 4: (Color online) (a) The logarithmic of second-order correla¬ 
tion functions (0) about the CCW-propagating intercavity field, 
versus DQE-field coupling strength parameter gj k, are calculated 
and shown with different mode-coupling strengths: black-solid line 
for J = 0; blue-dashed line for J — 10/c; red-dotted line for 
J = 20/^, and yellow-dot-dashed line for J = 40/^. (b) The second- 
order correlation function (0) of the CCW-propagating intercav¬ 
ity field is plotted as functions of both the mode-coupling strengths 
J//C and the DQE-field coupling strength parameter g/ i<i. The other 
system parameters for this simulation are: n = 407a, A = 0, 
£ = 7a, and 7a = 1 MHz. 


without the mode coupling in the WGMC are coupled to the 
DQE, the photon antibunching is observed around A = 0. 
Such antibunching in the weak DQE-field coupling regime 
results from an interference between the coherent light trans¬ 
mitted through the resonant cavities and the super-Poissonian 
light generated by photon-induced tunneling 11^ . Thus a bal¬ 
ance between the decay rate n of the cavity modes and the 
DQE-field coupling g is required. When the DQE-field cou¬ 
pling is not in the strong-coupling regime, the antibunching 
effect is weak, e.g., the minimum value of ^ccw(^) about 
0.47 at A = 0 for /i: = 2^, as shown in Eig. [3] However, 
when the mode coupling J is introduced, an obvious mini¬ 
mum and valley appear around A = 0, e.g., the minimum 
value ^ccw(^) reaches 0.006 at A = 0, in contrast to the 
case when there is no mode coupling. This strong antibunch¬ 
ing is due to the mode coupling, which will be discussed in 
the following paragraphs. 

To clearly show how the DQE-field coupling strength g af¬ 
fects the antibunching for different mode-coupling strengths 

J, including J = 0, 5'ccw(^) versus g is plotted in Eig.HJa) 
for different values of J. As shown in Eig.[4ja), for no mode¬ 
coupling, i.e., J = 0, we find: (i) a minimum value of 
^ccw(^) ~ around of g/i<i = 0.2, which is in the weak 
DQE-field coupling regime; (ii) the value of ^ccw(Q) 
creases with increasing g /k further; (iii) the value of (^) 
approaches one when g/ n> and the photon antibunching 

gradually disappears; (iv) the value of ^^rials to one 

and the photon antibunching does not occur when g/tv > 1, 
which is in the strong DQE-field coupling regime. This is 
counterintuitive with the fact that the better photon antibunch- 
ing (i.e., the smaller value of (0)) corresponds to a larger 

g. The physical mechanism for this effect has been explained 
in Ref. 11^ . Therefore, a balance between g and k is required, 
which might not be easy to achieve experimentally. Eor ex- 



EIG. 5: (Coloronline) Schematic diagram of the energy levels of 
the WGMC-DQE system corresponding to Fock states up to two- 
photons. The first photon and second photon indexes in the ket cor¬ 
respond to the photon number in the CCW propagating mode and the 
CW propagating mode, respectively. The index + (—) indicate the 
excited state (ground state) of the DQE. The long blue solid line at 
the bottom denotes the ground state of the whole system. The nar¬ 
rower short gray lines indicate single-photon states with decay rate 
tv, and the wider short gray lines at the top represent two-photon 
states with decay rate 2tv. Here we have neglected the decay of the 
DQE due to the fact that ^ 7a. The four green double-arrows 
represent the driving optical pump for the CCW propagating mode. 
The six blue double-arrows represent different energy level transi¬ 
tions induced by the DQE-field coupling. The four red duble arrows 
represent different energy level transitions induced by the CCW-CW 
mode coupling. 


ample, photon antibunching can occur only in a narrow re¬ 
gion around the dip, e.g., g/tz = 0.2. However, as shown in 
Fig-ffia), when mode-coupling is introduced, a better photon 
antibunching is achieved both in the weak and strong DQE- 
field coupling regimes. Eigure Ua) clearly shows that the 
better photon antibunching (smaller value of 5'ccw(^)) 
responds to the stronger mode-coupling strength in both the 

weak and strong DQE-field coupling regimes. The variations 

(2) 

of gccwi^) gj Kdio similar for different values of J. 

We now study how the statistical properties of the cav¬ 
ity fields are affected by the mode-coupling strength J and 
the DQE-field coupling strength g by plotting ^ccw(^) 
Eig. lUb). We find that a strong photon antibunching can be 
achieved both in the DQE-field weak-coupling and strong¬ 
coupling for the strong mode-coupling, e.g., J/ti > 5. Eor 
a given g, the larger J corresponds to the smaller value of 

5'ccw(^)’ ^ better photon antibunching is achieved. 

Similarly, for a nonzero mode-coupling strength J, the larger 

g corresponds to the smaller value of 5'ccw(^)’ ^ 

better photon antibunching is achieved. Erom Eig.[4l we also 
find that a nonzero mode-coupling strength J makes the bal¬ 
ance between g and ti unnecessary, to obtain a strong photon 
antibunching. This will be further discussed below. 
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C. Analytical solutions 

To better understand numerical calculations, we now study 
an analytical solution and compare it with numerical results. 
We assume that the driving field is very weak so that the total 
excitation number of the system is no more than two ll39l - l4^ . 
In this case, we use the ansatz 


\(p{t)) = Co,o,- |0,0, —) + If,0, —) + C'o,!,- |0,1, —) 

+ ^^0,0,+ |0,0, +) + C'2,0,- |2,0, —) + C'0,2,- |0,2, —) 

+ 1^5 ~) + ^1,0,+ |1: O5 +) + ^^0,1,+ |0,1, +), 

( 22 ) 


with Cij^± = Cij^±{t) , to calculate the steady state of the 
system. In the weak-driving limit, we have the relation 

<^ 0 , 0 ,- ^ C'o,!,-, Co, 0 ,+ 

^ C' 2 , 0 ,-,<^ 0 , 2 ,-, <^ 1 , 0 ,+: ^" 0 , 1 ,+. (23) 

Using the relation in Eq. ([23]) and combining with Eq. ([21]) . 
we have 

9giw(0) = (24) 

1^1,0,-1 

The coefficients (^ 2 , 0 ,- and Ui^o,- in the steady state can be 
obtained via the Schrodinger equation 

ih^=H\^), (25) 

and d\(p) /dt = 0. The effective non-hermitian Hamiltonian 
^ in Eq. ([25]) is given by 

= ft (a - (a^a + b^b) + hJ {a'<b + b^’a) (26) 

+ h az + h [g{a^ + b^)(J- + ea^ + H.c.] , 

in the zero temperature approximation. Here, the dephasing 
of the DQE has been neglected. 

Then, replacing H in Eq. ([25l) by the expression in Eq. ([26l) . 
we can obtain a set of linear equations, the detailed derivation 
can be found in Appendix A. By solving these linear equa¬ 
tions, we can obtain the coefficients, e.g. (^ 2 , 0 ,- and Ui,o,-, 

in Eq. ([22l) . then ^ccw(^) further calculated as 

2 1^2,0,-!^ _ l^ll^ 1 ^ 2 !^ 

with 

= ^p[^p^d + (AdAp — 2^^)(Ad + J)] + 

^2 = Ad(J + Ap) — 2^^, 

^3 = (^ + Ap)(J + Ap + Ad) — 2^^, 

A 4 = ApAd - . 


^ccw(^)) 


(27) 


(28) 

(29) 

(30) 

(31) 



FIG. 6: (Color online) The probability amplitude of the WGMC- 
DQE system with one-photon and two-photons versus the detun- 
ning parameter A//c are calculated and shown with different mode¬ 
coupling strengths: the black-solid curve for the one-photon proba¬ 
bility amplitude |(7i,o, -and blue-dashed curve for the two-photon 
probability amplitude |(72,o,-|^ with J = 0; the red-dotted curve 
for the one-photon probability amplitude |(7i,o,-|^ and the yellow 
dot dash curve for the two-photon probability amplitude |(72,o,-|^ 
with J = 20/5). The other system parameters for this simulation are: 
Hi = 407a, g = 207a, £ = 7a, and 7a = 1 MHz. 


Here, we define two detunnings by including the dissipation 
rates 


Ap = A-z|, (32) 

Ad = A-z^. (33) 

Eor comparison with the numerical calculations using the 
master equation, gcc-^^(O), calculated by analytical solution 
Eq. ([ 27 ]) . versus A/k, is also shown in Eig.[3] It is clear that 
the analytical solutions agree well with the numerical calcula¬ 
tions. 


D. Physical mechanism and comparisons between with and 
without mode coupling 

Using analytical solutions, let us further understand the 
physical mechanism of the mode-coupling-enhanced photon 
antibunching. In Eig. [5] we schematically show the energy 
level structure of the system up to two-particle excitations. In 
this case, nine states are involved, they are the ground state 
|0,0,0), three one-particle excitation states |1,0,—), |0,1,—), 
and |0, 0, -F), and five two-particle excitation states |2, 0, —), 
|0,2,—), ), |1,0,-|-), |0,1,-f). The connections be¬ 

tween different states, due to the driving field 5 and the cou¬ 
plings, characterized by J and g, are also schematically shown 
in Eig. [5] 
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We first consider the single-particle excitation. In this case, 
the CCW mode is pumped to the single-photon state and 
the system is in the state |1,0, —). When there is no mode 
coupling, i.e., J = 0, the probability amplitude 
which represents the probability to find the system in the state 
11, 0, —), is only determined by the decay rate n of the CCW 
mode and DQE-field coupling strength g. In this case and as 
shown in Fig. [5l there are two transition paths from the state 
|1,0, —) to the state |0,0, +) or the state |0,0, —). However, 
when the mode coupling is introduced by the scatter, | Ci,o, -1 ^ 
depends not only on n and g but also on the mode-coupling 
strength J. Thus an additional transition path between the 
state 1 1 , 0 , —) and the state |0, 1 , —) is introduced by the mode 
coupling, which may reduce the probability of the system at 
the state |1,0, —). 

When the CCW mode is further pumped to the two-photon 
state |2, 0, —). The additional transition path between the state 
|2, 0, —) and the state |1,1, —) is introduced due to the mode 
coupling, and this may also reduce the probability of the sys¬ 
tem in the state |2, 0, —). Therefore, when mode coupling J is 
introduced, both the single-photon probability |Ci,o,-|^ and 
two-photon probability |C 2 ,o,-|^ might be reduced, simulta¬ 
neously. Considering Eq. (EH), we now ask a natural ques¬ 
tion: which probability is larger, because only the ratio be¬ 
tween |C 2 ,o,-|^ and |Ci,o,-|^ makes sense. 

The probability amplitudes |Ci,o,-|^ and |C 2 ,o,-|^ versus 
detunning A/^c are plotted in Fig. 0 for J = 0 and J = 20/i:. 
We analyze Fig.[6|in two parameter regimes, i.e., |A|//^ < 0.6 
and |A|//i: > 0.6. In the first parameter regime and around 
A = 0, itis clear that both |Ci,o,-|^ and |C 2 ,o,-|^ are reduced 
when the mode coupling is introduced. They are changed as 
follows 

|C^2,o,-Ij=2o. = 10-®|C^2,o,-Ij=o> (34) 

and 

I^ia-IL20. = 10“"I^ia-Ij=o- (35) 

Considering the expression of 5(ccw(Q)’ 

^ccw(^) reduced by two orders of magnitude when the 
mode coupling is introduced, i.e., stronger photon antibunch¬ 
ing is achieved. When |A| /k > 0.6, both |C 2 ,o,-|^ and 
I Cl, 0 ,- 1 ^ are comparable for the case with mode coupling 
or without mode coupling, thus same for 

the cases with and without mode coupling in this parameter 
regime. 

We can also obtain the value of 5(ccw(^) analytical so¬ 
lution in Eq. (A9). When J = 0, we have 


1^1.0,-! 


2 _ 
J=0 ~ 


|ApAd - g^l 
|Ap|2 1-252+ ApAd|2' 


(36) 


Around A = 0, the ratio of |Ci,o,-|^ between J = 0 and 
J 7 ^ 0 is given 


1^1,0,- I J^Q 

1^1.0.-iLo 



(37) 



FIG. 7: (Color online) Logarithmic plot of second-order correla¬ 
tion functions CCW-propagating intercavity field as 

functions of both cavity field decay rate n and DQE-field coupling 
strength parameter g without optical mode coupling in (a) and with 
optical mode coupling rate J = 800 MHz in (b). The balance line 
between parameter g and /c is shown as a dash-dotted line in (a). 
The strong and weak DQE-field coupling regime are labeled up and 
below the dash line respectively. The other system parameters are: 
A = 0, £ = 7a, 7a = 1 MHz. 


Similarly, we can have 


l<^2,0 ,-lj^O 

|C^2,0,-|Lo 


4J6' 


Then, we can obtain 5ccw(0) 


as 


gccw(0)j?^o ^ -R 2 ^ 4^2 
5ccw(*))-t=0 -R2 J2 


(38) 


(39) 


When J = 20k, we find that 5ccw(0) ~ range 

near A = 0 , which is the same as numerical calculations. 

According to the above discussions, we conclude that the 
mode coupling results in the strong photon antibunching in 
the range near A = 0. To clearly show this, ^ccw(^) 
ther plotted as a function of k and g with J = 0 in Fig. [TJa) 
and J = 800 MHz in Fig. Hb). In contrast to the case of 
DQE-bimodal system without mode coupling the strong 
photon antibunching appears in the larger parameter range of 
the DQE-field coupling strength g when the mode coupling is 

introduced. Comparing this with no photon antibunching, i.e., 

(2) 

9ccwi^) ^ 1, in the strong DQE-field coupling regime when 
there is no the mode-coupling, we find that a strong photon 
antibunching, i.e., ^ccw(^) 10 “^, appears in the strong 

DQE-field coupling regime when the mode coupling is intro¬ 
duced. We also find that a balance between k, and g is nec¬ 
essary to obtain the stronger photon antibunching in the weak 
DQE-field coupling regime when J = 0. However, when the 
mode coupling J is introduced, such balance is not required. 


V. MODE COUPLING ENHANCED ROBUSTNESS TO 
PURE DEPHASING 

It is well known that dephasing with the rate 7 p cannot 
be neglected in some particular DQEs, e.g., the dephasing 
due to electron-phonon coupling play a very important role in 
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7p/7a 

FIG. 8 : (Color online) (a) The second-order correlation functions 
^ccw(O) CCW-propagating intercavity fields are plotted as 

functions of the cavity-field detuning A//c, with 7 p = 37a. The 
black-solid curve is for the cavity mode without coupling (J = 0 ) 
and red-dash-dotted curve is for strong mode-coupling (J = 10 /^). 
(b) The 5'ccw(0) versus the dephasing rate 7 p/ 7 a with different 
mode-coupling strength: J = 0 (black-solid curve), J = 6 (blue- 
dash curve), J — 10/^ (red-dotted curve) and J — 20a^ (yellow- 
dotted-dash curve) at A = 0. The other system parameters for these 
calculations are: n = 407a, g = 207a, £ = 7a, 7a = 1 MHz. 


quantum dot systems. Thus, we now study the pure dephas¬ 
ing effect and show how it affects the statistical properties of 
CCW propagating optical mode in the weak DQE-field cou¬ 
pling regime. 

We first study ^ function of the detuning A. 

Figure [3] has shown that the photon antibunching can be 
achieved around A = 0 for J = 0 when there is no dephasing. 
However, as shown in Fig.^a), when a pure-dephasing, e.g., 
7 p = 37a, is introduced, the photon antibunching effect dis¬ 
appears with all the values of ^q^(O) > 1 for J = 0. Thus, 
if there is no the mode coupling 11 ^ . the photon antibunching 
is not easy to be achieved when there is the dephasing. How¬ 
ever, when the mode coupling is introduced, e.g., J = 10/^, 
as shown in Fig. Oa), a strong photon antibunching can be 
achieved even when there is the dephasing. The minimum 
value of ^ccvv(^) corresponding to J = 10/^ is 0.06 at the 
point A = 0 in contrast to ^ccw(Q) = ^ J = 0. 

We further study how the dephasing affects ^ccw(^) 
weak DQE-field coupling regime in different mode-coupling 
strengths. Figure [ 8 tb) shows that increases with in¬ 

creasing 7 p and quickly approaches to one for J = 0. How¬ 
ever, when mode coupling is introduced (as shown in Fig.[ 8 tb) 
for three different mode-coupling strengths, e.g., J = Qn, 
J = 10K, and J = 20/^), the photon antibunching can still 
occur even with large dephasing. Interestingly, there is a min- 
imum value of ^ccw(^) versus of 7 p for each mode-coupling 
strength, and this minimum value represents the optimal case 
for photon antibunching. We find that photon antibunching 
can still be achieved when 7 p becomes very large. This means 
that the system is robust to dephasing when mode coupling 
is introduced. Figure [Sfb) also shows that larger value of the 
mode-coupling strength J provides a better photon antibunch¬ 
ing when the DQE has dephasing. 

In Fig.Oa), the ^ccw(^) plotted as a function of A and 
7 p with J = 0 and J 7 ^ 0, e.g., J = 20i<i. For J = 0, shown 



FIG. 9: (Color online) The second-order correlation functions 
^ccw(O) of the CCW-propagating intercavity fields are plotted in 
(a) as functions of both the cavity field detuning A//^ and DQE de¬ 
phasing rate 7 ^, without mode coupling J = 0 (the cyan-curved 
surface), and with mode coupling J = 20/c (the magenta-curved sur¬ 
face). (b) shows 5'ccw(0) versus both mode-coupling strength J/n 
and DQE dephasing rate 7 ^ at A = 0. The other system parameters 
are: /c = 407a, g — 207a, £ = 7 a, and 7 a = 1 MHz. 

as a cyan-curved surface, in Fig. [9ta), gccw(^) Quickly ap¬ 
proaches to one around A = 0 with increasing jp. This means 
that the antibunching is significantly affected by the dephas¬ 
ing. However, for J = 20k, shown as a magenta-curved sur¬ 
face, in Fig. [3a), ^ccw(ll) slowly increased with increas¬ 
ing 7 p around A = 0, and remains smaller than one up to 
Jp = 10 ja. Thus the system for the photon antibunching is 
very robust to the dephasing 7 p when the mode coupling is 
introduced. 

In Fig. [3b), ^ccw(ll) plotted as a function of J and 
7p. Figure [3b) shows that the larger mode-coupling strength 

J corresponds to the smaller value of ^ccw(ll) ^ given 
Jp, and thus a better photon antibunching. It is also obvious 

that all the values of gccw(^) smaller than 0.1 for all jp 
when J/k > 4. This reveals that strong photon antibunch¬ 
ing can occur even when a large dephasing jp is introduced 
(e.g., for 7p/7a = 20 and even larger). This is very differ¬ 
ent from the case without the mode coupling, i.e., the value 
of ^ccw(ll) approaches one even with a small dephasing, as 
shown in Fig. [3b) for J = 0. Therefore, we conclude that 
photon antibunching is robust to the phase dephasing of the 
DQE when the mode coupling is introduced. 

VI. ROBUSTNESS TO MODE MISMATCH AND 
COUPLING MISMATCH 

In previous sections, we mostly focused on the resonant 
interaction between the DQE and cavity modes and also the 
homogeneous couplings between two cavity modes and the 
DQE, i.e., uj = uja and ga = gb = g- Let us now analyze 
the effect of frequency mismatch and coupling mismatch on 
photon antibunching. 

We first study the effect of the mismatch between the fre¬ 
quencies UJ and cca of the cavity modes and the DQE on the an¬ 
tibunching when mode coupling is introduced. In Fig.[T0l for 
given = 0.5b:, the second-order correlation function 

^ccw (f^) plotted as a function of two detunings A = uj—uj^ 
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FIG. 10: (Color online) The second-order correlation functions 
^ccw(O) CCW-propagating intercavity field versus both the 

optical mode detuning A/k, and the DQE detuning Aa/n, for (a) 
J = 0 and (b) J = 20/^. The other system parameters are: 
= 407a, g — 207a, e = 7a, and 7a = 1 MHz. 


and Aa = cja — ujd, for J = 0 and also for J ^ 0, e.g., 
J = 2{)k. We find that ^ccw(^) reaches its minimum value 
when the driving field is resonant to the DQE, e.g., = 0. If 

we fix Aq = 0, the effect of the mismatch of the frequencies 
uj and cca on the antibunching can be clearly observed only 
by changing the detuning A. Compared with Fig. [I3a), in 
Fig.fTOlb). we find that increasing the mismatch | A|, the value 

of gccw(^) increases very slowly. The value of ^ccw(^) i^" 
creases to about 0.1 when A is changed from zero to a finite 
value for Aa = 0. So, when mode coupling is introduced, the 
photon antibunching can always be obtained in any detuning 
A when the driving field is resonantly applied to the DQE, 
i.e., when Aq = 0 . 

We then study the effect of the mismatch of the coupling 
strengths ga and gf^ on the antibunching when the mode cou¬ 
pling is introduced. In Fig. Ela) and Fig. fTlT b). for given 

^a = 0.5 k and A = Aa, we plot ^ccw(^) functions of 
A and g^ with J = 0 and J = 20 /^, respectively. Fig. [TlT a) 
shows that a small range around A = 0 and g\y = = 

O.S^a corresponds to ^ccw(^) ^ which means photon an¬ 
tibunching. However, in the other range, gQQ^^ (0) > 1, which 
corresponds to no photon antibunching. The minimum value 
appearing in Fig. EIu) for J = 0 is very sensitive to g\^ and 
A. We find that antibunching cannot be obtained even when 
the DQE-field coupling strength gij becomes large. Thus a 
balance between g\y and g^ is required to achieve photon anti¬ 
bunching when there is no mode coupling. 

When the mode coupling is introduced, e.g., J = 

Fig- nU b) shows that the value of decreases 

and then increases very slowly when increasing ^5 for A = 0 . 
The minimum value of 5'ccw(^) located at g\^ = k, which 
corresponds to the strongest antibunching. We find that in 
the weak DQE-field coupling regime, i.e., g\^ < k, the larger 
value of gi, corresponds to the smaller value of 5 'ccw(^) 
then better antibunching. However, in the strong-coupling 
regime, i.e., gi, > k, the value of ^ccw(^) increases with 
increasing gij, but the value of gccwi^) small, 

i.e., gccw(^) — 0-0005 in the given parameter regime. 

To further find the effect of the mismatch between g^ and g^ 







-4 


bJO -6 


1 

0.5 

0 0 g^/ 27 r 


FIG. 11: (Color online) The second-order correlation functions 
^ccw(O) of tbe CCW-propagating intercavity field are plotted as 
functions of both cavity field detuning A/k and DQE-field coupling 
rate gb/K (without mode coupling J = 0 and with mode coupling 
J — 20 /c) for p = 0.5 /^, as shown in (a) and (b), respectively. 
^QQw(O) is plotted as functions of coupling strength g^ and g^ (with¬ 
out mode coupling J = 0 and with mode coupling J = 20/^) for 
A = 0, as shown in (c) and (d), respectively. The other system pa¬ 
rameters are: k — 407a, £ = 7a, and 7a = 1 MHz. 
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on the photon antibunching, ^ccvs^(O) is plotted in Fig.fTTTc) 
and in Fig.[Tnd) as a function of Qq, and for A = = 0, 

with J = 0 and J 7 ^ 0, e.g., J = 20k,. For J = 0, Fig. fTTT c) 
clearly shows the antibunching can only occur in a narrow 
range which needs a balance between and as previous 
discussions. In most of the parameter range of ga and the 

value of ^ccw(f^) satisfies the condition gccw(^) — f ’ which 
means that there is no photon antibunching. However, when 
mode coupling is introduced, e.g., J = 20/i:, the balance be¬ 
tween ^fa and g\) is not required. Figure [TTId) also shows that 
the value of ^ccw(f^) increases when increasing both g^, and 
gi) for given other parameters. From Fig.fTTJd), we also find a 

very small value 10 “^ of gccw(^) achieved at g^j = k 

when ^a = 0, i-c., the CCW mode is decoupled from the DQE. 

Therefore, we conclude that the strong robustness gainst 
mismatches discussed above can be achieved when mode cou¬ 
pling is introduced. A balance, between ga and g^ for achiev¬ 
ing better photon antibunching for the case without mode 
coupling is not required when mode coupling is introduced. 
In particular, we find that the mismatch between ga and 
can play a positive role to produce photon antibunching when 
mode coupling is introduced. 


VII. CONCLUSIONS 
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Appendix A: Analytical solutions of the second-order 
correlation function with weak pump limit 

In the limit of a weak-driving field, the total excitation 
number of the system is assumed no more than two. Using 
the ansatz given in Eqs. (22) and (23) and combining with 
Eqs. (24)-(26), the coefficients in the steady state can be 
obtained via the Schrodinger equation by set 9 | (/?) = 0. 

With the zero temperature approximation and neglecting the 
pure dephasing of the DOE, we finally obtain the following 
set of linear equations 


In summary, we have studied photon antibunching in a sys¬ 
tem consisting of a bimodal WGMC and a DQE. Two modes 
in the WGMC are coupled via a scatterer. We mainly study 
the effect of the mode coupling on the photon antibunching, 
it is known that photon antibunching cannot occur in a sys¬ 
tem where a single-mode cavity field is weakly coupled to a 
DQE. However, we find that a very strong antibunching can be 
achieved in the weak DQE-field coupling regime, when mode 
coupling is introduced. 

Comparing with the case when the DQE is coupled to two 
modes of the bimodal cavity without mode coupling 
when mode coupling is introduced, we find that: (i) a strong 
photon antibunching can be achieved in a larger parameter 
regime including both the weak and strong DQE-field cou¬ 
pling regime; (ii) the minimum value of second-order corre¬ 
lation function can be reduced up to several orders of magni¬ 
tude even in the weak DQE-field coupling regime; that is, a 
strong photon antibunching can be achieved; (iii) the system 
is more robust to phase dephasing of the DQE; (iv) the sys¬ 
tem is also more robust to the frequency mismatch between 
the cavity modes and the DQE, as well as the mismatch of the 
coupling strengths between two modes of the cavity and the 
DQE. 

Our study shows that the system of two coupled cavity 
fields interacting with a DQE could be a good quantum de¬ 
vice for producing antibunched photons. Our studies support a 
platform to achieve single-photon sources on chip with robust¬ 
ness against fabrication imperfections and several mismatches 
in the system. 


0 

0 

0 

0 

0 

0 

0 

0 


+ ^(7o,o,+ + (Al) 


(A — i —)Co,i,- + JCi^o,- + 5'C'o,o,+? (A2) 

(A _ 2^)(7q gCi^Q_ -h (A3) 

(2A - m)C2,o,- + \/2 

-K\/2^Gi^o,+ + ? (A4) 

(2A - m)Co,2,- + V2+ y2^^Co,i,+(A5) 
(2A — m)(7i^i^_ -j- a/2JGo,2,- + V2JC2^o,- 
+^C'o,i,+ + ^C'i,o,+ + (A6) 

2A — i(—^" 1 , 0 ,+ + -^67o,i,+ 

+V2gC2,o,-gCi^i^- (A7) 

2A — i( — 2 ^^ ^^0,1,++ ^^^1,0,+ 

+ V^^Go,2,-- (A8) 


Due to the weak pump limit 5 ^ 0, we can assume Go,o,- ^ 
1, and one additional equation, namely eCi^o,- = 0, is ir¬ 
relevant to the problem. The Eqs. dATTi-dASt are now closed 
(i.e., eight equations for eight parameters). Thus, it is pos¬ 
sible to obtain the analytical solutions for all the coefficients 
Cij^±. However, the solution is cumbersome, so only Ci,o,- 
and 672,0,- are given 
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where 


Cl 

C 2 

C 3 


= -2g^ 





L = Ad + A2 + JA2 Ad 


<^ 1 , 0 ,- 

e(ApAd - 5 ^) 

(J - Ap)Ci ’ 

(A9) 

- 2 ApAd 5 (^ - 2JApg'^ + 5 ^, 

A A 

(A14) 

(A15) 


V2e'^L 

C 3 ’ 

(AlO) 

z 

Ad = A-iy. 

(A16) 


JAd 

\2 




(J + Ap) + Ad(J + Ap) — 2g , 

2(J-Ap)2CiC2, 


(All) 

(A12) 

(A13) 


We substitute Eqs. (IA9I) and (lAlOl) into Eq. (24), and then can 
derive the analytical solutions of the ^ccw(^) form of 
Eq. (27). 
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